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ABSTRACT 


In this study, the author develops a MATLAB simulation of 
area search with acoustic sensors modeled by the Poisson Scan 
model and the Lambda-Sigma (A—o) model. Detection time results 
are compared to those given by the much simpler Random Search 
formula. Random Search was found to closely approximate the more 


complex models if detection range was selected correctly. 





Guidelines for selecting the Random Search detection range were 


developed. 
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EXECUTIVE SUMMARY 


In this study, the author develops a MATLAB simulation of 
area search with acoustic sensors modeled by the Poisson Scan 
and the Lambda-Sigma (A—o) models. 

Both the Poisson Scan model and Lambda-Sigma (A-—o) model 


simulation results are found to be approximately exponentially 








distributed, which is consistent with the Random Search model. 


Thus, Random Search with the proper deterministic detection 


range R can be used to closely approximate simulation results 





obtained using the Poisson Scan and Lambda-Sigma (A—o) models. 

Both the Poisson Scan and Lambda-Sigma (A—oa@) detection 
time results vary with A, o, and R(50). The author constructs a 
regression model based on the simulation results and finds an 


approximate linear relationship between R(50) and the best-fit 





R for reasonable A and o. The gradient of the regression line 
depends on the values of A and o. Thus, it is possible to 


estimate the best-fit R from problem parameters A, a, and R(50) 





and then use Random Search, instead of Monte Carlo simulation, 





to predict th ffectiveness of area search. 
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I. INTRODUCTION 


A. BACKGROUND 


Random Search is one of the most well-known and used models 





for area search. However, because of its simplicity, it also 








has inherent limitations. For example, it assumes that the 


searcher and target search areas are identical and that the 








searcher uses a perfect cookie-cutter sensor. Captain Gi Young 
Kim!, ROK Air Force, studied these issues previously and 


suggested generalizations to the Random Search formula for 





cases where the searcher and target areas are not identical. 








However, there are remaining issues to address, particularly 


how well the Random Search model can approximate probabilistic 





sensors, such as those modeled by the Poisson Scan and 


Lambda-Sigma (A-—o) models. 


B. RESEARCH QUESTIONS 


The primary research questions are the following: 


1. When using the Poisson Scan model and the Lambda-Sigma 
(A-—o) model, will detection times be approximately 
exponentially distributed, as predicted by the Random 
Search model? 











2.Can a non-cookie cutter sensor be replaced with a 
cookie-cutter sensor with the same sweep width and 
maintain approximately the same search performance? 


3. Are there significant differences in the detection results 
generated by simulations using the Poisson Scan and 
Lambda-Sigma (A—o) detection models? 





1 Gi Young Kim, “Development and Testing of a New Area Search Model with 
Partially Overlapping Target and Searcher Patrol Areas,” Master’s thesis, 
Naval Postgraduate School, Monterey, California, 2008. 


il 


Cc. THESIS ORGANIZATION 


The thesis is organized as follows: Chapter II reviews the 
Poisson Scan and Lambda-Sigma (A—@) acoustic detection models; 


Chapters III and IV present a MATLAB simulations of these two 








models; Chapter V develops the best-fit R, which is the 
equivalent cookie-cutter detection range for the Poisson Scan 


model and Lambda-Sigma (A—o) model; Chapter VI summarizes all 





results and recommends future studies. 


II. OVERVIEW OF ACOUSTIC DETECTION MODELS 


A. BACKGROUND 





In this chapter, we review the Passive Sonar Equation (PSE) 
and use it to develop two models for passive acoustic detection 
— the Poisson Scan model and the Lambda-Sigma (A—@a@) model. 


B. THE PASSIVE SONAR EQUATION 


The passive sonar equation? is an audit of energy flow 





between a source and a sonar receiver. All units are in 
decibels (dB). 
The Passive Sonar Equation(PSE) is 


Ls — Ny — DNL — DT =SE, 





where 
Ls; = signal level of source, 
Ny = propagation loss between the source and 
receiver, 


DNL = detection noise level at the receiver, 





DT = detection threshold = signal level at receiver 





required to achieve a Py of 0.5, 


SE = signal excess. 





The left side of Passive Sonar Equation is the signal 





available at the receiver, and right side is the signal required 


for a probability of detection(P,) of 0.5. 





2 The references for this section are the unpublished lecture notes of 
Professor James N. Eagle, “Acoustic Detection Models,” 2009. 
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b, Figure of Merit(FOM) and R(50) 


Figure of merit (FOM) and R(50) are used as measures of 


ffectiveness(MOE) for sonar detection performance. FOM is 





the maximum propagation loss a sonar systemcan absorb in dB and 
still produce a 0.5 probability of detection. R(50) is the 


maximum range between the source and receiver resulting in a 





0.5 probability of detection. 


FOM = L;— DNL — DT 


= the transmission loss between source and receiver 





resulting in Py=0.5. 


R(50) = the maximum range where L; =FOM. 


2. Acoustic Modeling Assumptions 


e Mean SE=SE=FOM-—Ny. That is, the Passive Sonar 
Equation is used to calculate mean signal excess at the 
receiver. 

e SE~N(SE,o7), where o is typically 3-9dB. We assume the 
Signal excess is a normally distributed random variable 
with the mean given by the Passive Sonar Equation, and 
0 is specified. 


e Detection occurs if and only if SE2>0dB. This is a 





threshold crossing model. Detection occurs when the 
random SE is nonnegative and at no other times. And 


probability of detection (called here instantaneous 











probability detection, IPD) is P(SE20O). 








SE ~ N(SE, 0°) 
IPD = P(SE > 0) = ®(SE/c) 


Figure 1. IPD when SE <0 


3. Implications 


IPD(SE) = P(SE =>0), where SE~N(SE,o7) 
= @(SE/o). 


where ® is the cumulative normal distribution function. This 





relationship between SE and the probability of detection is 





illustrated in Figure 1. 


Cc. THE FIXED SCAN AND POISSON SCAN MODELS 
1. Fixed Scan Model 


The Fixed Scan model? assumes that probabilistically 





independent detection opportunities occur at regularly spaced 


time intervals. It is attractive analytically, since there is 








a simple formula for calculating the probability of detection 











over any time interval [0, t]. Let S be the set of detection 





opportunity times T within the interval, and let 


IPD(t) = ®(SE(t)/o). 


3 Alan R. Washburn, Search and Detection, 4th ed. 3-4p. 
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be the detection probability for a scan containing time T, where 


the signal excess SE(t) is acknowledged to depend on time. Then 








the cumulative detection probability over the interval is 


PUT St) = F(t) = 1—[pes(1 — IPD()). 


In other words, there will be a detection at [0, t] unless 





detection fails at every one of the independent scans within 





the interval. However, the Fixed Scan model suffers from its 
dependence on an arbitrarily selected origin of time. The 
Poisson Scan model was developed, in part, to address this 


issue. 


2. Poisson Scan Assumptions 





e Independent detection opportunities occur at Poisson 
times with rate A(units : 1/time). 

@e Probability of detection at time t, given that a 
detection opportunity exists at t, is IPD(t) = ®(SE(t)/o). 

3:5 Development 


From the Poisson Scan assumption, the instantaneous 


detection rate at time t is 


y(t) = AB(SE(t)/o). 


Then it follows from the properties of the non-homogeneous 


Poisson Process that 


P(T <t) = Fy(t) =1—exp (-Af_, @(SE(t)/o)ds). 


In this model, 1/A is the mean time between independent 
detection opportunities (glimpses). Probabilistic independence 
is assumed for each glimpse, with A specifying the average 


glimpse rate. By varying A, the sonar can be made as effective 








(A larger) or ineffective (A smaller) as desired, subject to 


sensor or operational constraints. The sonar system can only 








call detection at the discrete glimpse times, and otherwise is 


assumed to be processing previously received data. 


D. THE LAMBDA-SIGMA (4-0) MODEL 


In the Poisson Scan model, detection can occur at time 
t only if a detection opportunity occurs at t, and an independent 
draw from SE~N(SE, 07) is greater than or equal to 0. 

A potentially more realistic model would assume that SE(t) 


is a random, continuous function of time t (that is, a 





continuous - time stochastic process), and that detection 


occurs at any time t where SE(t)20. 


1. (A—o@) Model Assumptions and Development 


Let 
SE(t) = SE(t) — X(t), 


where SE(t) is the mean signal excess at time t, which is 


computed(as in the Poisson Scan model) from the PSE as 


SE(t) = FOM(t) — Ny (t). 


X(t) is a 0-mean stochastic process. More specifically, X(t) is 


a (A-—o) jump process, described below. 


Each random sample path of X(t) is a step function, where 


the duration of each step (the time between jumps) is 





exponentially distributed with mean 1/4, and the height of each 


step (the value of X(t)) is normally distributed with mean 0 and 


variance oa”. 









initial 
detection 


X ~ N(0,07) 





time. r 





SE(t) loss of contact 
< Er > 
T ~ exp(A) 
Figure 2. (A—o) Detection Process Example 


As illustrated in Figure 2, detection occurs whenever 
SE(t) => X(t). X(t) can be thought of as the time-varying sum of all 


random components of FOM and Ny. 


2. Computing Py when SE is Constant 


Assume SE(t) =SE =k, t € [0, tT]. We will develop an expression 
for the probability of initial detection occurring between time 


0 and t, which we call F,(t). 


F,(t) = P(X(t) < k, forsomet € [0, tT]) 
=1-P(X(t)>k, forallt € [0, t]) 
=1- j=) P(X() > k, forall t € [0, t]|i jumps occur in [0, T]) 
x P(i jumps occur in [0, T]) 
- 44 (Atye _ 
= 1-Yo(1 — P(k/o))#1 Pe 


((1-9(k/o)]at)! 
i! 


=1-[1- ®(k/o)Je™ Yo 
=1-[1- O(k/o)Je Mel PHO, since YRo— = e* 
=1-[1- O(k/a)Je*Mer/ole 


=1-[1- BGSE/a)je MK | 


F(t) 





[Peso sie 0) = (En 








= 
time t 
tp 
= 
time t 
Figure 3. (A—oa) PDF and CDF Functions for Time of 


Initial Detection T when SE is Constant 


At time 0, the probability of detection is ®(SE/o). 








The mean time to detection, E(T), and the density function, 


fr), for time of initial detection are 


1 — @(SE/o) 
A®(SE/a) ” 


fr(t) = A®(SE/o)[1 — (SE /a)Je4PS4/9t, for t>0. 
Also, P(t =0)= @(SE/a). 


E(T) = 


Note that this probability distribution has a discrete portion 





(T=0) with units of probability (unitless) and a continuous 


portion (T>0) with units of probability/time (or, 1/time). 
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III. SIMULATION OF POISSON SCAN MODEL 


A. DESCRIPTION OF POISSON SCAN SIMULATION MODEL 
a I Characteristics of the Searcher 


The searcher is not allowed to search for a target outside 





the search area. In addition, the searcher’s initial position 


is uniformly distributed inside of the search area. After that, 





the searcher selects his course randomly, independent of the 








target’s movement. The course change event is determined by 
Poisson process with rate A,. Independent scan opportunities 
also occur according to a Poisson process, but with rate d;. Also, 
it is assumed that signal transmission loss4# follows either the 
spherical spreading law (N,=20logigr) or the cylindrical 


spreading law (N, =10logi;9r). The searcher speed was fixed at 





15 nm/hour . 


2. Characteristics of the Target 


Like the searcher, the target is not allowed to move 








outside of the search area and has an initial position uniformly 











distributed over the search area. The logic of the target 











movement is the same as that of the searcher; that is, the target 
has its own, independent Poisson process with course change rate 


A,. The target speed was fixed at 5nm/hour. 


4 Robert J. Urick, Principles of Underwater Sound, 3rd ed. 101p. 
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B. COMPUTER ALGORITHM 


5 


Br 


Input [Units] 


Number of simulation replications, Nreps = 500. 





Maximum simulation time, tmg, = 150[hour]. 
The length of search area in X direction, l, = 150[nm]. 
The length of search area in ¥Y direction, l, =150[nm]. 


Searcher speed, V [nm/hour]. 





Target speed, U [nm/hour]. 

Searcher’s scan rate, A, [Glimpses/hour]. 
Searcher’s course change rate, A, [1/hour]. 
Target’s course change rate, A; [1/hour]. 
The unit time of simulation, At [hour]. 

The size of search area, As=k, Xl,. 
Figure of Merit, FOM [dB]. 

Transmission loss, MN, = 20log(Distance) [dB]. 
Signal Excess, SE=FOM—N,,. 


Variance of Signal Excess, o [dB]. 


Functioning of the Poisson Scan Simulation 


When a new replication begins, the initial positions of 








the searcher and the target are chosen from a uniform 


distribution over the search area. The initial headings are also 





drawn uniformly between 0 and 2m radians. The subsequent course 


changes for the searcher and target occur according to Poisson 


processes with rates A, and A,, respectively. According to 


12 





Captain Gi Young Kim’s thesis®, the recommended value for the 


course change rate is A=2V/,/l,xl,, which implies that on 











average, two course change events occur during the time required 





for the searcher to go from edge to edge in the search area. 








When the searcher or target encounters an area boundary, 
a random reflection occurs. After a reflection, the new course 
in radians is Uniform_Random(C, —0.5,C,+0.5), where C, is the 


perpendicular course from the reflection boundary. This scheme 





was recommended and tested by Captain Gi Young Kim® to achieve 
an approximate uniform distribution of the searcher and target 
tracks in the search area. 


In the Poisson Scan model, there are two ways to simulate 











detection times. One uses the instantaneous probability 











detection at specific time t, which follows the Poisson process 





with rate A. This is because in the Poisson Scan model, detection 


opportunities occur according to a Poisson process with rate 








A. So, the mean time between independent detection opportunities 





is 1/A. The distance between the searcher and the target is 





calculated at specific time t. Then, the mean signal excess(SE) 
is determined by using the Passive Sonar Equation(PSE) allowing 
the calculation of IPD(SE). The next step is generating a 
standard uniform random number and comparing it with the IPD(SE). 
If Uniform Random(0,1) < IPD(SE), then a detection occurs at that 


specific time t. 





2 Gi Young Kim, “Development and Testing of a New Area Search Model with 
Partially Overlapping Target and Searcher Patrol Areas,” Master’s thesis, 
Naval Postgraduate School, Monterey, California, 2008. 


6 Gi Young Kim, KOREA AIR FORCE CAPTAIN, graduated from the Naval 
Postgraduate School in 2008. 
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Mean Signal Excess[dB] 
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Figure 4. Simulation when calculating distance and SE at 
specific time t 





The other way to simulate detections is to check for a 


detection at each At. The distance between the searcher and the 








target is calculated at every time step At, and the mean signal 


excess is computed. 
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Figure 5. Simulation when calculating distance and SE at 
every At 


After that, we compute the instantaneous probability of 
detection IPD(SE) = @(SE(t)/o) and detection rate y(t) =A®(SE(t)/c). 
At each time step At, if Uniform Random(0,1)<y(t)At, then a 


detection occurs. 
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Time 


Figure 6. Fr(t) for different simulations 


As illustrated in the simulation results of Figure 6, the 
results of these two simulation methods are very similar, and 
theoretically should be identical. 

For this thesis we used the At simulation method. In order 
to approximate a continuous simulation, At should be small. 
However, too small a At results in too many calculations. In the 


next section, we answer the question, “How small should At be 
to produce accurate results without creating excessively long 


simulation runs?” 
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CDP 





Time 


Figure 7. F7(t) for various values of At 


3. Output 


R(50) 


Extensive simulation experimentation showed that At = ae 








was the largest At resulting in reliable and repeatable results. 


Thus we move the searcher half the distance of the R(50) at each 





time step. Typical results are shown in Figure 7. 


In order to find best-fit R, which is the equivalent range 


of a cookie-cutter sensor to the Poisson Scan Model, the author 





experimented with various values of dA, o, and FOM. 
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FOM in Poisson Scan model 
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Fr(t) for various Sigma in Poisson Scan model 
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Figure 10. 
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F,(t) for various Lambda in Poisson Scan model 





The Green line represents the Poisson Scan model formula. 


P(T <t) =F,(t) =1-— exp caf @(SE(s)/o)ds). 


s=0 


The red line represents the formula resulting from a random 
search with a cookie-cutter sensor that has best-fit R range. 
P(T <t) = F(t) = 1— exp (—2RVt/A,), 


~~ In(4 7 Fr(tmax)) * A, 
2V tas 


R= 


J 





where V is the approximate mean relative speed between target 


and searcher determined as follows: 





= 1 7 
v=-| V2+ U2 —2VU cos@ dé. 
T Jo 








The blue line represents the results of the simulation. 











The results show that for the model parameters examined 
here, the Poisson Scan model gives detection times that are 


approximately exponentially distributed, as does the Random 








Search model. The conclusion is that, the Random Search model 


using the proper detection range R can closely approximate 


detection results given by the Poisson Scan model. 
A potentially more accurate procedure for determining R 
would be to plot —In(1—F;(t)) versus t and to solve for the 


best-fit slope y. Then, 
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IV. SIMULATION OF LAMBDA-SIGMA (A—o) MODEL 


A. DESCRIPTION OF LAMBDA-SIGMA (A—o) SIMULATION MODEL 


The Lambda-Sigma (A—a) detection model’ is distinguished 





from the Poisson Scan model by allowing detection at any time 
during the scan, rather than only at specific scan times. This 
model assume that SE(t) is a random, continuous function of time 


t. Then, 


SE(t) = SE(t) — X(t), 


where X(t) is a 0-mean stochastic process where the duration of 


each step is exponentially distributed with mean 1/A, and the 











height of each step is normally distributed with mean 0 and 


variance o%. Detection occurs at any time t where SE(t)>0. 





7 The references for this section are the unpublished lecture notes of 
Professor James N. Eagle, “Acoustic Detection Models,” 2009. 
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B. COMPUTER ALGORITHM 
1. Input [Units] 


e Number of simulation replications, Nrep, = 500. 


e Maximum simulation time, tm, = 150[hour]. 





e The length of search area in X direction, l, =150[nm]. 











e The length of search area in Y direction, lL, = 150[nm]. 





e Searcher speed, V=15[nm/hour]. 

e Target speed, U=5[nm/hour]. 

e Searcher’s course change rate, A, = 2V/ Jl, x ly [1/hour]. 
e Target’s course change rate, A, = 2U/,fl, X ly [1/hour]. 
e The unit time of simulation, At = R(50)/2 x V [hour]. 

e The size of Search area, A,=l, Xl,. 

e Figure of Merit, FOM [dB]. 

e Transmission loss, Ny = 20log(Distance) [dB]. 


e Signal Excess, SE =FOM—N,,. 





e Jump rate of stochastic process, A, [times/hour]. 
e Standard deviation of jump in stochastic process, o[dB]. 
2% Functioning of the Program 


With the Lambda-Sigma (A-—o) model, the duration of each 


step is exponentially distributed with mean 1/d,, and the height 








of each step is normally distributed with mean 0 and variance 


ao. 
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X(H[dB] 





Time 


Figure 11. Example of (A—o) model simulation 


In Figure 11, the blue line represents one realization of the 





O-mean stochastic process X(t), and the red line represents mean 


signal excess. At each time step At, if SE(t)>X(t), then a 
detection occurs. In the Figure 11, the first detection occurs 


at 55 hours. 


3. Output 


In order to find the best-fit R, which is the equivalent 


range of a cookie-cutter sensor to that of the Lambda-Sigma 





(A—o@) model, the author experimented with various values of 


model parameters A, o, and FOM. 
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Figure 12. F;(t) 


for various FOM in 


(A—oa) model 


24 



























Sigma = 0dB Sigma = 1dB 
og cl 
06 eo 
5 as 05 
ae a 
ae a : Random Search formula 
Be A ——: Simulation result 
5 @0 700 760 5 20 700 150 
= Time = Time 
R = 4.61 y = 0.0063 R=5.01nm y = 0.0069 
Sigma = 2dB Sigma = 3dB 
og O9-F 

















oO 50 100 150 



































Time 2 aie 
R =5.11nm y = 0.0070 R =5.89nm y = 0.0081 
Sigma = 4dB Sigma = 5dB 
i 1 
og 09 
% osl & os 
O4r o4 
a A 0.2 
95 Ei 700 760 95 Ea 700 760 
= Time = Time 
R =6.56nm y = 0.0091 R =7.84nm y = 0.0107 
Sigma = 6dB Sigma = 7dB 
= os 
fs 04 
°5 =o 700 760 700 150 
Time cliente 


R= 10.08nm y = 0.0138 





R =12.68nm y = 0.0174 








Figure 13. 


F7(t) for various Sigma in 


(A—@) model 
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Figure 14. Fr(t) for va 


rious Lambda in (A—o) model 
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The red lines in Figures 12-14 show Random Search results 








with a cookie-cutter sensor having best-fit R range. The blu 








line represents the results of the simulation. As was the case 


with the Poisson Scan model, simulation results using the 





Lambda-Sigma (A—o) model showed an exponential time to initial 











detection. Thus, Random Search with the correct deterministic 








detection range R canbe used to closely approximate simulation 








results obtained using the stochastic Lambda-Sigma (A—o) model. 


The task that remains is to estimate from problem parameters 


the appropriate value of R. 
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V. ESTIMATION OF BEST-FIT R 


A. LATERAL RANGE CURVE 


ole Description of the Lateral Range Curve 


In order to estimate the best-fit R, the author simulated 
lateral range curves for sensors with the same parameters as 
the Poisson Scan model and Lambda-Sigma (A—oa@) model. Sweep 
width W is defined to be the area underneath the lateral range 


curve. 
W= | l(x)dx, 


where l(x) is the lateral range curve. In the cookie-cutter 


sensor with detection range R, 


(x) = { 1,2 € [-R, R] 


0, otherwise 


therefore, R=W/2. 





-L/2 -Rn searcher Rm Li2 


Figure 15. Sweep Width Interpretation 
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2. Computer Algorithm of Lateral Range Curve 


In the lateral range curve simulation geometry, shown in 





Figure 15, target lateral ranges are uniformly distributed over 
-L/2 to L/2 and target tracks are straight. Except for this 
characteristic, input elements in each lateral range curve 
Simulation are the same for both the Poisson Scan model and 
the Lambda-Sigma (A-—o) model. The lateral range curve is 
determined by the ratio of the number of detections that occur 
to number of replications at each lateral range from the 


searcher. 


Probability of detection 
° ° 
= 


° 
a 
T 


01 











0 1 1 1 L 
0 5 10 6 vt} 3 Ki} =.) “a 5 a 
Range 


Figure 16. Simulation result of lateral range curve 


Then, the area underneath the lateral range curve is sweep width, 


which can be computed numerically. 
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3. Output of Lateral Range Curve 


The following tables compare the half-sweep width computed 





from lateral range curve simulations to the best-fit R obtained 





from Poisson Scan and Lambda-Sigma (A-—o) simulations, for 


various FOM, dA, and o values and assuming spherical spreading. 























FOM{[dB] 74 77 80 84 87 90 
Rennes 1.40 | 2.74 | 4.94 | 8.24 | 13.08 | 18.15 
model [nm] 
(A—«a) model[nm] 2.80 3.97 Set 8.81 12.81 19.07 
Lateral range half 
sweep width of P-S 2.14 3.15 6.31 10.26 16.15 28.42 





model [nm] 
Lateral range half 
sweep width of 3.07 4.49 6.65 114331 16.87 25.63 
(A—o) model [nm] 















































Table 1. Comparison of R between Poisson Scan model, (A-o) 
model and half sweep width of each model when 
o =3dB, A=2——,V=15—™ pas ™ 
hour hour 








hour 





Sigma[dB]| 0 1 2 3 4 5 6 7 
Poisson Scan 
model [nm] 
(A-o) model[nm] 4.61 | 5.01 | 5.11 | 5.89 | 6.56 7.84 10.08 | 12.68 
Lateral range half 
sweep width of P-S | 3.68 | 4.29 | 5.12 | 6.31 | 8.07 | 10.62 | 14.56 | 20.78 
model [nm] 
Lateral range half 
sweep width of 9430) | 5.550) 5:2.98>| 26.64 | 77270 9.36 11.86 | 15.42 
(A-—o) model [nm] 






























































Table 2. Comparison of R between Poisson Scan model, (A-o) 
model and half sweep width of each model when 
FOM = 804B, A=2—_,V=15™ ,y=5™ 


hour’ hour 
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A[1/hour] 0.5 ih ie 2 245 3 
Poisson Scan 
model [nm] 
(A—«o@) model[nm] 5.03 5.58 5.86 5.89 5.99 6.06 
Lateral range half 
sweep width of P-S| 0.89 2.98 4.89 6.32 7.40 8.23 
model [nm] 
Lateral range half 
sweep width of oe 13 6.07 6.36 6.64 6.84 7.02 
(A—o) model[nm] 












































Table 3. Comparison of R between Poisson Scan model, (A-o) 
model and half sweep width of each model when 
FOM = 80dB, o = 3dB, V=15—",U=5-— 











hour’ hour 


These results showed that: 





e The best fit R is strongly dependent on FOM, A, and o for 
both the Poisson Scan and Lambda-Sigma (A—oa) models. 


e The Lambda-Sigma (A—o) best fit R somewhat exceeds the 
Poisson Scan best fit R, and this might be due to the 
Lambda-Sigma (A-—o) model starting with a positive 
probability of detection at time 0. 





e The best fit R generally increases with increases in FOM, 
A, and o for both the Poisson Scan and Lambda-Sigma (A-—o) 
models 








e The lateral range curve model produces half-sweep width 
values which can significantly exceed the best fit R 
values produces by the Poisson Scan and Lambda-Sigma (A-—o) 
models 














Tt is not clear why the lateral range procedure produced 
detection range estimates not matching well with those of the 
Poisson Scan and Lambda-Sigma (A—o) models; but it is possible 


that the lateral range assumptions (infinite, straight-line 





paths with uniformly distributed closest points of approach) 


were not well enough met in the area search simulation. 
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B. (50) 


5 Description of R(50) 





Another way to estimate the best-fit R is to use R(50), 





which is the maximum range between the searcher and target 





resulting in an instantaneous detection probability of 0.5. In 
other words, R(50) is the maximum range where mean signal excess 
SE is zero. For this analysis, the author assumes that signal 
transmission loss follows a mixture of spherical spreading 
(Ny =20logigr) and cylindrical spreading (Ny, = 10log,)7r) 

Specifically, we use Ny, =15logigr. The author also assumes that 


the source level of the target Ls; is in the range 84 to 95dB, 





detection noise level at the receiver DNL is 45dB, and 











detection threshold DT is -15dB. Therefore, by the passive 


sonar equation, FOM is in the range 54 to 65dB. 


2. Compute R(50) 


With the author’s assumptions, the FOM and R(50) can be 





calculated. 


SE = FOM —N,, 
= FOM —15logigr. 


Setting SE=0dB, Pz=0.5, and 1nm=1852m, we determine that 





FOM 
10 15° 


R(50) = Fa52 ° 


This value can be compared to the best-fit R values obtained 


from the Poisson Scan and Lambda-Sigma (A—o) simulations. 
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a. Output of R(50) 


For acoustic detection models, a reasonable 4A is 1 to 2hr7} 





and a reasonable o is 2 to 4dB. The author used regression to 


examine how well the computed R(50) values estimated the 





best-fit R values obtained from the Poisson Scan and 


Lambda-Sigma (A—o) simulations. 

















T r T L T : T r T ui T 


6 
R(50)[nm] 


Figure 17. Regression result of Poisson Scan model 
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4- 
07 T : T r T r T r T r T 
2 4 6 8 10 12 
F(50) [nm] 
Figure 18. Regression result of Lambda-Sigma (A—o) model 


As shown in Figures 17 and 18, the relationship between 
R(50) and best-fit R is approximately linear. In addition, the 
slope of each regression line strongly depends on the values 


of A, o, and the model being used. 





























o A 1 [hr *| 2 [hr] 
2dB R = 0.85 x R(50) — 1.71 R = 1.10 x R(50) — 1.67 
3dB R = 1.06 x R(50) — 2.17 R = 1.44 x R(50) — 2.41 
4dB R = 1.49 x R(50) — 3.05 R = 2.13 x R(50) — 3.79 
Table 4. Regression result of Poisson Scan model at each A and 
oO 
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oe A Lhe | 2[hr-*] 
2dB R = 1.05 x R(50) — 0.16 R = 1.14 x R(50) — 0.23 
3dB R = 1.26 x R(50) — 0.44 R = 1.50 x R(50) — 0.93 
4dB R = 1.65 x R(50) — 1.09 R = 2.15 x R(50) — 2.27 
Table. 5. Regression result of Lambda-Sigma (A-—a) model at each 
A and a 


As illustrated in Tables 4 and 5, it is possible to estimate 
the best-fit R from problem parameters R(50), A, and o. Then 


this best-fit R can be used with the Random Search formula below 








to further estimate the area search probability of detection 


by time ¢. 


2RVt 
P(T <t)=F,;(t)=1-e 4 
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VI. CONCLUSIONS 


A. CONCLUSIONS AND RECOMMENDATIONS 


There are two primary contributions of this thesis. The 





first is the demonstration that initial detection times for area 
search simulations using both the Poisson Scan and Lambda-Sigma 
(A-o) acoustic detection models are approximately 


exponentially distributed, allowing the simulation results to 





be closely approximated by the venerable Random Search formula. 





And the second contribution is the observation that the best-fit 
cookie-cutter detection range used in the Random Search formula 
can be accurately predicted using the simulation model 
parameters A, o, and the R(50) detection range. 

In this thesis, it is assumed that acoustic signal 
transmission loss follows either spherical spreading or 
cylindrical spreading. A potentially more realistic model could 


be developed by using actual propagation loss data. 
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APPENDIX 


A. MATLAB CODE 


as 


Poisson Scan Model (I) 


Nreps=500; Snumber of simulation replications 
tmax=150; Smax. simulation time (hr) 


1x=150; 
ly=150; 


search area length in x direction (nm) 
search area length in y direction (nm) 





V=15; Ssearcher speed (nm/hr) 

U=5; Starget speed (nm/hr) 

lams=(2*V)/sqrt(lx*ly); tsearcher course change rate (1/hr) 
lamt=(2*U) /sqrt (1x*ly); Ssearcher course change rate (1/hr) 


laml=2; 
sig=3; 
FOM=60; 


os 
(o) 


searcher looking rate (1/hr) 
signal excess variance (dB) 
figure of merit (dB) 





R50=10% (FOM/15) /1852; %R(50) (nm) 
dt=R50/ (2*V); Sdelta t (hours) 
Xs=zeros(1,tmax/dt+1l); %initialize x-position to zero(searcher) 





Ys=Xs; 
Cs=Xs; 
Xt=Xs; 
Yt=Xs; 
Ct=Xs; 
T=0:dt 
A=1x*ly 
Cc 

f 


os 
(o) 
o 
io) 
x 
io) 
os 
) 


x 
(o) 


, 





initialize y-position to zero(searcher) 
initialize searcher course to zero 
initialize x-position to zero(target) 
initialize y-position to zero(target) 
initialize target course to zero 





stmax; tsimulation time vector 


Ssearch area 


umDet=zeros(1, tmax/dtt+1); %initialize cumulative detection state 


or n=1:Nreps %main simulation loop 
xs=rand*1lx; %initial searcher and target x and y positions 
ys=rand*ly; 
xt=rand*1x; 





ll ct 





DGAKMAKKXATtTAAK 





, 








t=rand*ly; 

s=rand*2*pi; Sinitial searcher course 
=rand*2*pi; Sinitial target course 
O; 


Sset simulation time to 0 


dex=1; %initialize time index to l 





index)=xs; %save initial searcher x position 
index)=ys; %save initial searcher y position 


) 
) 
index)=cs; %Ssave initial searcher course 
) 
) 


index)=xt; tsave initial target x position 
index)=yt; %save initial target y position 





index)=ct; %save initial target course 
ance=zeros(1, tmax/dtt+1l); tSinitialize distance between target and 








=zeros(1, tmax/dtt+1); %initialize signal excess 
Gamma=zeros(1, tmax/dtt+1); %initialize instantaneous probability of 





detection 
Detection=zeros(1, tmax/dtt+1l); % initialize detection vector 
for t=l:tmax/dt Sinner loop 
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tindex = tindex+1l; tupdate simulation time index 
if rand<lams*dt; cs=rand*2*pi; end 

if Xs(tindex-1)<0; cs=(rand-0.5); end 

if Xs (tindex-1)>lx; cs=pi+(rand-0.5); end 

if Ys(tindex-1)<0; cs=pi/2+(rand-0.5); end 

if Ys(tindex-1)>ly; cs=-pi/2+(rand-0.5); end 
if rand<lamt*dt; ct=rand*2*pi; end 

if Xt (tindex-1)<0; ct=(rand-0.5); end 

if Xt (tindex-1)>(1x); ct=pit+(rand-0.5); end 

if Yt (tindex-1)<0; ct=pi/2+(rand-0.5); end 

if Yt (tindex-1)>(ly); ct=-pi/2+(rand-0.5); end 





Xs (tindex) = Xs(tindex-1)+V*dt*cos(cs); SUpdate x and y positions 
Ys(tindex) = Ys(tindex-1)+V*dt*sin(cs); 

Cs (tindex)=cs; 

Xt (tindex) = Xt (tindex-1)+U*dt*cos(ct); SUpdate x and y positions 
Yt (tindex) = Yt(tindex-1)+U*dt*sin (ct); 

Ct (tindex) =ct; 


Distance (tindex-1)=sqrt ( (Xs (tindex-1) -Xt (tindex-1)) .*%2+ (Ys (tindex-1)-Yt 
(tindex-1)).%*2); 
SE (tindex-1) =FOM-15*1log10 (Distance (tindex-1) *1852) ; 
Gamma (tindex-1)=laml*normcdf (SE (tindex-1)/sig) ; 
if rand <= Gamma (tindex-1) *dt; 
Detection (t: (tmax/dt+1))=1; 








end 
end Sinner loop (time increasing from 0 to tmax) 
CumDet = CumDet + Detection; 
end Souter loop (simulation replications) 
Probability=CumDet/Nreps; 
plot (T,Probability,'b-'), axis([0,tmax,0,1]) 
xlabel('Time', 'Fontsize', 12), ylabel('CDP', 'FontSize', 12) 
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2. Poisson Scan Model (II) 


Nreps=500; Snumber of simulation replications 
tmax=150; %Smax. simulation time (hr) 

1x=150; Ssearch area length in x direction (nm) 
ly=150; Ssearch area length in y direction (nm) 
V=15; Ssearcher speed (nm/hr) 


U=5; 


Starget speed (nm/hr) 


lams=(2*V)/sqrt(lx*ly); tsearcher course change rate (1/hr) 
lamt=(2*U) /sqrt (1x*ly); Ssearcher course change rate (1/hr) 
laml=2; %tsearcher looking rate (1/hr) 


sig= 
FOM= 
R50= 








3; SSignal excess variance (dB) 
60; Sfigure of merit (dB) 
10% (FOM/15) /1852; SR(50) (nm) 


dt=R50/(2*V); delta t (hours) 

Xs=zeros(1,tmax/dt+1); Sinitialize x-position to zero(searcher) 
Ys=Xs; Ginitialize y-position to zero(searcher) 

Cs=Xs; Sinitialize searcher course to zero 


ll cot ct 


QO: 





OF 


searcher 


Y 
G 
T 
A=1lx*ly; sSsearch area 
Cc 
ig 





Xt=Xs; Sinitialize x-position to zero(target) 
=Xs; Sinitialize y-position to zero(target) 
t=Xs; Ginitialize target course to zero 





dt:tmax; simulation time vector 





umDet=zeros(1, tmax/dtt+1l); %initialize cumulative detection state 


n=1:Nreps %Smain simulation loop 

xs=rand*1lx; %initial searcher and target x and y positions 
ys=rand*ly; 

xt=rand*1x; 

yt=rand*ly; 

cs=rand*2*pi; Sinitial searcher course 

c 








t=rand*2*pi; Sinitial target course 
t=0; tset simulation time to 0 
tindex=1; %initialize time index to 1 




















Xs (tindex)=xs; save initial searcher x position 

Ys(tindex)=ys; ssave initial searcher y position 

Cs(tindex)=cs; Ssave initial searcher course 

Xt (tindex)=xt; %Ssave initial target x position 

Yt (tindex)=yt; sSsave initial target y position 

Ct (tindex)=ct; Ssave initial target course 

Distance=zeros(1, tmax/dt+1l); Sinitialize distance between target and 








SE=zeros(1, tmax/dtt+1l); initialize signal excess 
Gamma=zeros(1, tmax/dtt+1); %initialize instantaneous probability of 








detection 


Detection=zeros(1, tmax/dtt+1l); % initialize detection vector 
for t=l:tmax/dt Sinner loop 

tindex = tindex+l; tupdate simulation time index 

if rand<lams*dt; cs=rand*2*pi; end 

if Xs (tindex-1)<0; cs=(rand-0.5); end 

if Xs (tindex-1)>1lx; cs=pit+(rand-0.5); end 

if Ys(tindex-1)<0; cs=pi/2+(rand-0.5); end 
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Ys(tindex-1)>ly; cs=-pi/2+(rand-0.5) 
rand<lamt*dt; ct=rand*2*pi; end 

Xt (tindex-1)<0; ct=(rand-0.5); end 
Xt (tindex-1)>(1x); ct=pit(rand-0.5); 
Yt (tindex-1)<0; ct=pi/2+(rand-0.5); 


if Yt(tindex-1)>(ly); ct=-pi/2+(rand-0. 
Xs (tindex) = Xs (tindex-1)+V*dt*cos (cs); 
Ys(tindex) = Ys(tindex-1)+V*dt*sin(cs); 
Cs (tindex)=cs; 

Xt (tindex) = Xt (tindex-1)+U*dt*cos (ct); 
Yt (tindex) = Yt(tindex-1)+U*dt*sin(ct); 
Ct (tindex) =ct; 

if rand<laml*dt 








7 end 


end 

end 

5); end 

SUpdate x and y positions 


sUpdate x and y positions 


Distance (tindex-1)=sqrt ( (Xs (tindex-1) -Xt (tindex-1)) .*%2+ (Ys (tindex-1)-Yt 


(tindex-1)).%*2); 
SE (tindex-1) =FOM-20*1log10 
Gamma (tindex-1) =normcdf (SI 
if rand <= Gamma (tindex-1); 
Detection (t: (tmax/dt+1))=1; 





D 
( 


= 
G 
py 








end 
end 
end Sinner loop (time increasing from 0 to 
CumDet = CumDet + Detection; 


end Souter loop (simulation replications) 
Probability=CumDet/Nreps; 
plot (T,Probability,'b-'), 
xlabel('Time', 'Fontsize', 


axis([0,tmax,0,1]) 
12), ylabel('CDP', ' 
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istance (tindex-1)*1852); 
tindex-1)/sig); 


tmax) 


FontSize', 12) 


ce Lambda-Sigma (A—0) Model 


Nreps=500; Snumber of simulation replications 
tmax=150; max. simulation time (hr) 
1x=150; Ssearch area length in x direction (nm) 
ly=150; Ssearch area length in y direction (nm) 
V=15; Ssearcher speed (nm/hr) 
U=5; Starget speed (nm/hr) 
lams=(2*V)/sqrt(lx*ly); tsearcher course change rate (1/hr) 
lamt=(2*U) /sqrt (1x*ly); Ssearcher course change rate (1/hr) 
lamda=2; Sduration of each step (1/hr) 
sig=3; Sheight of each step (dB) 
FOM=60; % figure of merit (dB) 
R50=10% (FOM/15)/1852; %SR(50) (nm) 
dt=R50/ (2*V); Sdelta t (hours) 
Xs=zeros(1,tmax/dt+1); %initialize x-position to zero(searcher) 
Ys=Xs; Sinitialize y-position to zero(searcher) 
Cs=Xs; Sinitialize searcher course to zero 
Xt=Xs; Ginitialize x-position to zero(target) 
=Xs; Sinitialize y-position to zero(target) 
=Xs; Sinitialize target course to zero 
=Xs; Sinitialize step function 
O:dt:tmax; Ssimulation time vector 
=lx*ly; Ssearch area 
umDet=zeros(1, tmax/dt+1); tinitialize cumulative detection state 
or n=1:Nreps %main simulation loop 
xs=rand*1lx; Sinitial searcher and target x and y positions 
ys=rand*ly; 
xt=rand*1x; 
yt=rand*ly; 
cs=rand*2*pi; Sinitial searcher course 
ct=rand*2*pi; Sinitial target course 
sp=randn*sig; %Sinitial step 
t=0; tset simulation time to 0 
tindex=1; %initialize time index to 1 














a €F UT 














Y 
Cc 
S 
T 
A 
Cc 
E 





























Xs (tindex)=xs; save initial searcher x position 
Ys (tindex)=ys; save initial searcher y position 
ti =Ccs; % initial 

Cs (tindex)=cs; %tsave initial searcher course 

Xt (tindex)=xt; %tsave initial target x position 

Yt (tindex)=yt; Ssave initial target y position 

Ct (tindex)=ct; Ssave initial target course 

SP(tindex)=sp; Ssave initial step 

Distance=zeros(1, tmax/dt+1); %initialize distance between target and 
searcher 

SE=zeros(1, tmax/dt+l); initialize signal excess 


Detection=zeros(1, tmax/dtt+1l); %initialize detection vector 








for t=l:tmax/dt Sinner loop 
tindex = tindex+1l; %tupdate simulation time index 


SP (tindex) =SP (tindex-1); 
if rand<lamda*dt; SP(tindex)=randn*sig; end 
if rand<lams*dt; cs=rand*2*pi; end 
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if Xs(tindex-1)<0; cs=(rand-0.5); end 

if Xs (tindex-1)>1lx; cs=pit+(rand-0.5); end 

if Ys(tindex-1)<0; cs=pi/2+(rand-0.5); end 
if Ys(tindex-1)>ly; cs=-pi/2+(rand-0.5); end 





if rand<lamt*dt; ct=rand*2*pi; end 

if Xt (tindex-1)<0; ct=(rand-0.5); end 

if Xt (tindex-1)>(1x); ct=pit+(rand-0.5); end 

if Yt(tindex-1)<0O; ct=pi/2+(rand-0.5); end 

if Yt (tindex-1)>(ly); ct=-pi/2+(rand-0.5); end 








Xs (tindex) = Xs(tindex-1)+V*dt*cos(cs); SUpdate x and y positions 
Ys(tindex) = Ys(tindex-1)+V*dt*sin(cs) ; 

Cs (tindex) =cs; 

Xt (tindex) = Xt (tindex-1)+U*dt*cos(ct); sUpdate x and y positions 
Yt (tindex) = Yt(tindex-1)+U*dt*sin(ct); 

Ct (tindex) =ct; 


Distance (tindex-1)=sqrt ( (Xs (tindex-1) -Xt (tindex-1)) .*%2+ (Ys (tindex-1)-Yt 
(tindex-1)).%*2); 
SE (tindex-1) =FOM-15*1log10 (Distance (tindex-1) *1852) ; 
if SP(tindex-1) <= SE(tindex-1); 
Detection (t: (tmax/dt+1))=1; 








end 
end Sinner loop (time increasing from 0 to tmax) 
CumDet = CumDet + Detection; 
end Souter loop (simulation replications) 
Probability=CumDet/Nreps; 
plot (T,Probability,'b-'), axis([0,tmax,0,1]) 
xlabel('Time', 'Fontsize', 12), ylabel('CDP', 'FontSize', 12) 
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